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Cumulation of High-current Electron Beams: Theory and Experiment 

S.V. AnishchenkoQ V.G. Baryshevsky0 N.A. Belous, A.A. GurinovichU 
E.A. Gurinovich, E.A. GurnevichU and P.V. Molchano\0 
Research Institute for Nuclear Problems 
Bobruiskaya str., 11, 220030, Minsk, Belarus. 

A drastic cumulation of current density caused by electrostatic repulsion in relativistic vacuum 
diodes with ring-type cathodes is described theoretically and confirmed experimentally. The dis¬ 
tinctive feature of the suggested cumulation mechanism over the conventional one, which relies on 
focusing a high-current beam by its own magnetic field, is a very low energy spread of electrons in 
the region of maximal current density that stems from a laminar flow profile of the charged-particle 
beam. 

PACS numbers: 84.70.-|-p, 52.59.Mv 


I. INTRODUCTION 

The pioneer research into high-current electron beams 
dates back to the 30ies of the past century [l|. For the 
lack of equipment and tools affording the generation of 
high-power charged-particle beams under terrestrial con¬ 
ditions, the researchers mainly focused their attention on 
theoretical consideration of astrophysical problems 0 . 

The first high-current electron beams with the power 
from several of gigawatts to several of terawatts ob¬ 
tained three decades later made a revolution in the cumu¬ 
lation research. This became possible primarily through 
two remarkable achievements in experimental physics: 
First, Dyke and colleagues experimentally obtained cur¬ 
rent densities as high as 10® amperes/cm^ from the mi¬ 
croprotrusions of the metal cathode placed in a strong 
electric field; second, the dielectric breakdown data re¬ 
ported by J. Martin and colleagues provided the 

potential for developing high-voltage pulse generators. 

Self-focusing of higlucurrent electron beams with their 
own magnetic fields [lll[l4| provided the charged-particle 
beam intensities as high as ^ ITW/cm^, thus enabling 
the laboratory investigation of the extreme state of mat¬ 
ter. The expectation was that by cumulation of high- 
current beams, the deuterium-tritium targets would be 
compressed and heated to ignition so as to initiate ther¬ 
monuclear reactions and thus accomplish pellet fusion 

mm- 

Though the initially set goal of developing a pellet fu¬ 
sion was not achieved, high-current electron beams found 
successful applications in other fields of physics JTM3. 
They are used for research in radiation phy sics [18|. gen¬ 
eration of high-power microwaves [l^ , collective ac¬ 
celeration of ions [2l|, , and pumping gas lasers . 

Nonlinear phenomena originating from the high-current- 
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beam interaction with self- and external electromagnetic 
fields figure prominently in all these processes. 

This paper considers one of such phenomena, which is, 
in fact, an alternative mechanism of high-power electron 
beam cumulation. This mechanism occurs in relativistic 
vacuum diodes with a ring-type cathode. Even though 
this phenomenon has been experimentally observed for 
years, it still lacks a consistent explanation. Our main 
task here is to provide a theoretical description of this 
cumulation mechanism and the experimental verification 
thereof. We will show in the subsequent sections that the 
explosive electron emission changing the emitting surface 
of a high-current diode is paramount for the cumulation 
process. 

This paper is arranged as follows: We shall first give 
quite a detailed description of the phenomenon of ex¬ 
plosive electron emission; then we shall describe the cu¬ 
mulation mechanism of the electron beam that was re¬ 
vealed during modeling the relativistic vacuum diode op¬ 
eration with the self-developed computer code (see 
the underlying computer algorithm is given in the Ap¬ 
pendix). In conclusion, we shall report the experimental 
results that confirm the described cumulation mechanism 
for high-current electron beams. 

II. ELECTRON EMISSION 

High-current electron beams are generated in relativis¬ 
tic vacuum diodes, composed of a cathode and an an¬ 
ode, through explosive electron emission (see fig.l). The 
mechanism of the explosive electron emission is as fol¬ 
lows [2B - [^ : once the voltage is applied across the elec¬ 
trodes of the relativistic diode, the field-emission current 
is emitted from the cathode surface [2^-[3l| , which is the 
electrons tunneling from metal into a vacuum under the 
influence of the electric field. The electrons moving in the 
metal heat the cathode surface. The microscopic electric 
field near the cathode is nonuniform because of the sur¬ 
face defects of the conductor. Particularly, the field near 
the microprotrusion tips is sufficiently greater than the 
average one, which causes rapid heating of the tips that 
explodes as the specific energy density rises to ~ 10^ J/g. 
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When the macrofield becomes as large as 1 MV/cm, 
the time delay td of the explosion is as small as 1 ns. 
Each microexplosion is accompanied by thermionic emis¬ 
sion from the surface of the cathode flare — conducting 
plasma expanding at a speed v ^ lO'^m/s (see Fig.2). 
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FIG. 1: Relativistic vacuum diode with a ring-type cathode 



FIG. 2: Explosive electron emission 

It is noteworthy that td is particular sensitive to the 
condition of the cathode surface. It was shown in [s^ 
that dielectric inclusions on the cathode surface result 
in the excessive increase of the field-emission current. 
B.M.Cox and W.T.Williams [s^ reported experimen¬ 
tally observed high electric field strengths near dielec¬ 
tric inclusions on the cathode surface with the local field 
strength in the vicinity of the inclusions being several 
hundreds times greater than the average field strength 
in the cathode-anode gap. The dielectric inclusions, as 
well as overall surface defects, naturally leads to the time 
spread td in the cathode flare formation in different mi¬ 
croregions of the cathode. 

Dielectric inclusions not only initiate explosive electron 
emission, but also sustain it |34| . The matter is that the 
activity of each emitting center is accompanied by the 
ion flow to the cathode (Fig.2). Dielectric inclusions in 
the vicinity of the emitting center are charged by the ion 


current resulting in the breakdown and the formation of 
new emitting centers. Another mechanism of cathode 
plasma formation is associated with the field desorption 
of the atoms absorbed on the cathode surface [s^. This 
occurs in the regions where the local electric field exceeds 
~ lO^V/cm. Collisional ionization of the desorbed atoms 
from the field-emission current results in the formation 
of plasma layer on the cathode surface. 

Expansion of the conducting plasma of the cathode 
flares leads to screening of the nearby regions on the 
cathode surface by a strong electric field. The analy¬ 
sis given in shows that the characteristic radius 

of the screened region is 

Tscr- 5- 10^ ( 1 ) 

where U[Y] is the applied voltage, /[A] is the cathode 
flare current, and /iac[cm] is the cathode-anode gap. For 
U = 400 kV, i = 10 A, and hac = 1-5 cm, the charac¬ 
teristic radius of the screened region, according to (P), is 
Tscr = 0.15cm. 

Thus, the number of the explosive-emission centers Ng 
occurring simultaneously on the surface of the cathode 
of radius Tc = 3 cm can be estimated at Ne = r’^jr'^cr = 
60. The number of concurrent explosive-emission centers, 
Nf., indicates the degree of inhomogeneity of the beam’s 
transverse structure. 


III. CUMULATION MECHANISM 

As it has been stated in the previous section, explo¬ 
sive electron emission begins with the formation and ex¬ 
pansion of cathode flares. Explosive electron emission is 
most intense from the cathode protrusions, particularly 
from the cathode’s inner edge (Fig.3). Coulomb repul¬ 
sion causes the charged particles to rush to the region 
free from the beam. As a result, the accelerated motion 
of electrons toward the anode comes alongside the radial 
motion to the cathode’s symmetry axis. As a result, the 
high-current beam density increases multifold on the axis 
of the relativistic vacuum diode as compared to the av¬ 
erage current density in the cathode-anode gap. The re¬ 
ported cumulation mechanism are described for the ring- 
type cathode (the circular cathode with a hole, which 
coincides with the cathode exis). It should be noted that 
the cumulation mechanism doesn’t depend on the hole 
position, i.e. whether or not the hole axis coincides with 
the cathode axis. 

Figure 4 shows the results of simulations: the dose 
absorbed by the anode. The assumed parameters of the 
cathode were as follows: cathode radius 3.0 cm, cathode- 
anode gap 2.0 cm, and the radius of the inner hole 0.8 
cm. The maximum value of the accelerating voltage pulse 
was taken equal to 360 kV and its duration - to 330 
ns. The simulated current density in the region of the 
central spot on the anode at the moment corresponding 
to the maximum accelerating voltage was as large as 1.0 
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U=360 kV, r)j=0.8 cm, ftac=2.0 cm, n^=\A cm ^ 



Cathode 

flares 


FIG. 3: Cumulation caused by electrostatic repulsion 

kA/cm^, being 5 times greater than the average current 
density of the high-current diode. Thus, the simulation 
result indicate the electron-beam cumulation on the axis 
of a high-current diode with a ring-type cathode. 

The undeniable advantage of this cumulation mecha¬ 
nism over a conventional one based on the self-focusing 
of a high-current beam by its own magnetic field is a 
very low energy spread of particles in the region of the 
maximum current density due to the laminar flow of 
charged particles (Figure 5). In contrast, under the con¬ 
ditions of self-focusing of the beam by its own magnetic 
field, the flow current becomes appreciably turbulent, 
and the charged particles acquire a significant momen¬ 
tum spread (39j. The electron flow in this case is like a 
compressed relativistic gas with electron temperature of 
the order of the voltage applied across the diode (l^ . 


IV. EXPERIMENTAL RESULTS 
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FIG. 4: Absorbed dose (simulation results) 


U=360 kV, r>j=0.8 cm, cm, n^-\A cm ^ 



2 4 6 S 10 12 14 16 


X, cm 


To investigate the cathodes and obtain the information 
about electron beam parameters, we used a nanosecond 
pulse-periodic electron accelerator with a compact SF6- 
insulated high-voltage generator (HVG) as a power sup¬ 
ply providing pulsed voltage up to 400 kV in ~ 30 Ohm 
resistive load with half-height duration of 130 ns and rise 
time of 30 ns [40l |. 

To obtain integrated full-sized imprints of the electron 
beam, we used a radiochromic dosimetry film (techni¬ 
cal specification TU 2379-026-13271746-2006) (U placed 
3 mm behind the anode mesh made of stainless steel 
(the geometrical transparency of the mesh was 0.77); the 
cathode-anode gap was 2 cm. The dosimetry film en¬ 
abled us to obtain information about the total absorbed 
dose due to the passage of charged particles [4l |. 

After the exposure, the transmission scanning of the 
dosimetry film using the optical filter was made with EP¬ 
SON Perfection VlOO Photo scanner. The distribution of 
the absorbed dose (and hence the beam’s energy density) 


FIG. 5: Electron flow (simulation results) 


over the beam cross section was derived unambiguously 
from the scanned images and the dose-response calibra¬ 
tion curve. 

Our first experiments showed that the flow of charged 
particles on the axis was so intense that it burned the 
film through (see Fig.6). For this reason, we placed 70 
/rm-thick aluminium foils in front of the dosimetry film 
to decrease the absorbed dose (see Figs. 7-10). 

This enabled us to cut off the flows of both the weakly- 
relativistic electrons produced at the voltage pulse decay 
and the cathode plasma. The experiments conducted 
with one, two, and three foils demonstrated that a sharp 
increase in the absorbed dose remains in the center. This 
means that the particle flow consists at the beam axis of 
high-energy electrons. In the experiments with three foil 
layers cutting off all electrons whose energy is less than 
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FIG. 6: Electron beam imprint without foils 
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FIG. 7: Experimental scheme 


250 keV, the absorbed dose in the center was almost four 
times as large as the average dose across the beam cross 
section, showing a good agreement with the simulation 
results. 

Let us note here that both the simulation and the ex¬ 
periments were performed at maximum accelerating volt¬ 
age ^ 400 kV and the cathode-anode gap equal to 2 cm. 
The estimates show that with the voltage increased to 
2 MV and the cathode-anode gap decreased by a factor 
of 5 it is possible to achieve the beam intensity of the 
order of ITW/cm^ required, for example, to study the 
extreme states of matter and to do research into inertial 
confinement fusion. 


V. CONCLUSION 

In this paper we described the cumulation mechanism 
of a high-current beam in a relativistic vacuum diode 
with a ring-type cathode. The basis of this cumulation 
mechanism is electrostatic repulsion of electrons from the 


FIG. 8: Electron beam imprint with one foil 
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FIG. 9: Electron beam imprint with two foils 

explosive-emission plasma on the inner edge of the cath¬ 
ode. The simulated values of current density and beam 
intensity equal to 0.36GW/cm^ and 1 kA/cm^, respec¬ 
tively, qualitatively agree with the experimental data. 

A very low particle energy spread in the region of 
maximum current density that is due to laminar flow of 
charged particle is the distinctive feature of the described 
cumulation mechanism over a conventional one relying on 
focusing the high-current beam by its own magnetic field. 

Appendix A: Simulation of high-current beams 

In simulating of the electron beam dynamics under 
the conditions of nonuniform explosive electron emis¬ 
sion, it is necessary to consider the expansion of the 
cathode plasma emitted from separate explosive emis¬ 
sion centers. Self-consistent simulation of particle mo- 
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FIG. 10: Electron beam imprint with three foils 


tion in self- and external electric and magnetic fields is 
usually performed using the particle-in-cell method 
Si in a quasi-stationary approximation [d^, [d^ . Quasi¬ 
stationary approximation applies when the field param¬ 
eters in high-current diodes change slowly, and the dis¬ 
placement currents and induction fields are neglected. 

The simulation is performed in the Cartesian, as well 
as in cylindrical coordinates. The electric fields and par¬ 
ticle motion are computed in the Cartesian coordinate 
system, while the magnetic fields - in cylindrical. Spatial 
dimensions of the cells on the mesh for the field calcula¬ 
tions are set equal, i.e., AX = AT = AZ = AR. The 
system is assumed to be axially symmetric, but the nu¬ 
merical simulation of particle motion is performed in 3 
dimensions, which is of principal necessity for a proper 
consideration of the electron emission nonuniformity. 


1. Particle-in-cell method 


The particle-in-cell method has been developed for 
simulating the multiple phenomena in different fields of 
physics: vacuum electronics [48l - [5ll| . plasma physics 


^|, hydrodynamics [62|, magnetic hydrodynamics 1631 . 
64\, astrophysics [h^, and semiconductor physics [^.IhTj. 

This method consists in the representation of real flows 
of charged particles (electrons, protons, and ions) as a set 
of macroparticles, each containing a large number of real 
charge carriers. Every macroparticle, normally located in 
a single cell, has a certain attributed spatial distribution 
of mass and charge. Depending on the charge and the 
location of the macroparticle, a certain contribution to 
the charge and current densities is attributed to its near¬ 
est nodes on the spatial mesh using the weighting proce¬ 
dure. Using a similar procedure, one can find the forces 
acting on the macroparticle, knowing the magnitudes of 
the fields in the nodes located in the close proximity to 
the particle. Substituting the value for the force into 
the finite-difference analogues of relativistic equations, 
we may find new locations and momenta of macropar¬ 
ticles. For a full description of the system of fields and 
particles, these procedures are completed with charged 
particle injection into and removal from the computa¬ 
tional region. 

It is the numerically realized injection of charged par¬ 
ticles emitted from the surface of the expanding cathode 
flares that constitutes the novel and important feature of 
the developed code. 

A typical program cycle based on the particle-in-cell 
method consists of six operations: computation of coor¬ 
dinates and momenta of particles, injection and removal 
of particles, computation of the current and charge densi¬ 
ties, and computation of the electric and magnetic fields. 
Now let us proceed to a detailed description of each pro¬ 
cedure realized in our code. 



FIG. 11: Particle-in-cell 


2. Electric fields 

The electric field strength in the quasi-stationary ap¬ 
proximation in the Coulomb gauge is related to the scalar 
potential (j) governed by the Poisson equation: 

A(j) = —inp, (A-1) 

E = (A2) 

The analysis given in showed that the Jacobi iterative 
method fits best to solve the Poisson equation for plasma 
dynamical problems. It is well known that the iter¬ 
ative methods requiring that the approximate value of 
the potential at the first iteration step be specified are 
slowly converging, because the initial distribution usually 
deviate appreciably from the exact solution of the finite- 
difference analogue of the Poisson equation. In plasma 
dynamical problems the situation is basically different, 
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which stems from a more appropriate selection of the 
initial approximation at the first iteration step: the grid 
potential magnitudes obtained at the previous time step 
are taken as the first approximation [^. As a result, 
the entire iteration process at each time step reduces to 
one-three iterations, requiring much less time than, say, 
the computation of new coordinates and positions of the 
particles. 

Thus, with the Jacobi iterative method the final- 
difference analogue of the Poisson equation has the form: 


After we find the potential, the electric field is found im¬ 
mediately from [toI] 

/ J-P- AP AP AP 

pn { ~ 9i,j+i,k ~ 9i,j-i,k 

i,j,k - -1 2AX ’ 2Ay 

'Pi,j,k + l ~ 9i,j,k-l \ 

2AZ )' 


3. Magnetic fields 


in,s 

^i+l,j,k 


-2C 






jn,s 

^i,j + l,k 


+' 


^i,j,k+l 


AA2 

-2C^r+' 


Ay2 


■ 71,S 

'i,j,k- 


AZ^ 


= —4np. 


'i,j,k ! 


ri,j,k 


4P~^ 
ri^j^k 5 


■ 71,S 

^i,j — l,k 


(A3) 


where s is the iteration number and n is the time step number. 
The iteration in (IA3II occurs until the residual 
becomes less than The value of which is the 

norm of the (/>"•" matrix is computed by formula 


In considering motion of relativistic charged particles, it 
is fundamentally important to take account of the self- and 
external magnetic-field effect on the electron-beam dynamics 
in a high-current diode. As we are concerned with axially 
symmetric configurations, we shall proceed to cylindrical co¬ 
ordinates. By virtue of axial symmetry, we shall assume that 
the magnetic field H and the current density j are indepen¬ 
dent of the azimuth angle 6. 

In the absence of the axial field Hp, Hg becomes the only 
magnetic field component and is related to the current density 
jz and the current I running through the cathode by the 
Stokes theorem 




(A4) 


The parameter e in fact determines the error of the finite- 
difference Poisson equation solution. 

To solve the Poisson equation, we need to complete it with 
the boundary conditions. We used the Dirichlet boundary 
condition implying the specified potentials on the cathode 
{(f> = Uc) and the anode {4> = Ua)- At the edge of the com¬ 
putational region in the cathode-anode gap we took the log- 
arithmical distribution of the potential [7ll| 


(l}=Uc + 


{Ua - Uc) Injr/rA 
ln{ra/rc) 


(A5) 


that exactly describes the change in (j) in the gap between two 
infinite cylinders. Obviously, the distribution of the potential 
in a high-current diode will approach (IA7II if the boundary of 
the computational region stated here is located at a consid¬ 
erable distance from the electron-emitting surface. We shall 
find the grid density Pij^k with the linear weighting procedure 
attributing a certain contribution to p^p^k coming from eight 
nodes nearest to the particle located at {xa,ya, Za) 


^Ptlk = qP - wPl - 0(1 - wZ)/AV 
Ap“+i,yfc = q.wTil - 0(1 - O/AP 
ApOi./c = qp - wrO(l - w^)/AV 
PTlk+i = qp - wpi - wJ)wl/AV 
Ap“0+i.fc = g.«'>7(l-0/AP ^ ^ 

ApOi,.+i = qp- w-)wpi/AV 

Ap“+A,,. = Q„u>r(l - 0«^fc/AP 

A cx.n n n n i \ -t t 

Api + Ij + I_fc + 1 = qaWi WjWk/AV. 

Here AP = AXAYAZ w" = (^(a;" - Xij,k)/AX, (p" - 

yiPlAYAz^^-ZiPIAz). 


He 


4n{I{r,z) + 2n jz(ri, z)ridri) 
2-Kcr 


(AS) 


Here the current I{r, z) contains the contributions coming 
from all electrons injected into the points with coordinates 
> 2. 


4. Cathode plasma expansion 

The center of explosive emission is formed on the cathode 
when the electric field strength exceeds a certain threshold 
value Ear that depends on the condition of the electrode sur¬ 
face. Then the cathode flare begins to expand at a speed 
Veee ~ 2 • 10® cm/s in many directions. Without deliber¬ 
ate control over the surface microstructure, the emission cen¬ 
ters are chaotically located about the cathode surface. The 
mean distance between them is determined by the size of the 
screened area dacr = 2rsc7- ~ 3 mm (see O), knowing which 
we can easily estimate the characteristic density rieee of the 
explosive-emission centers rieee = \/'Kr\ar ~ 14 cm“^. 

For the purposes of high-current diodes simulation, we shall 
assume that the emission regions are formed in the cathode 
nodes where the electric field exceeds Ear with the probability 
neea{AXAYAZ)^^^. When the electric field becomes greater 
than Ear, the cathode flare expands at a constant speed Vees 
in every direction from the emission region. Each cathode 
flare is the source of electrons. The active cathode flare emits 
thermionic emission current which is many times as large as 
the current limited by the beam space charge, and so we can 
speak about the unlimited thermionic emission resulting in 
practically zero field on the surface of the expanding cathode 
plasma. The assumption about the cathode screening enables 
us to appreciably simplify the numerical computation of the 
charged-particle kinetics, sparing ourselves the need to simu¬ 
late fast processes just in the emission region. Such simulation 
would require a high space-time resolution due to the small¬ 
ness of the Debye length and high values of plasma oscillation 
frequency [fi^. 
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At each time step, we inject the charged particles into 
plasma-occupied nodes that have in the vicinity at least one 
node free from the conducting material. The magnitude of 
the injected charge Qij,k is found from the relation 

( rpn _ rpn rpn _ rpn 

■^xi-\-l,j,k -^xi — l,j,k ■^xi,j-\-l,k -^xi,j — l,k 

2AX 2Ay 

^ 2I^Z ) dvr 

(A9) 


Let us note that the charge is injected if Qi,j,k < 0. 


5. Motion of charged particles 


Numerical integration of relativistic equations of motion is 
the most time-consuming procedure of all, that is why it is 
paid special attention to by the developers of the codes. In 
a nonrelativistic case, the most widely used is the leapfrog 
scheme [H, HJ] 


+ F2/ST, 


^+1/2 

-«n+l -fn I Pot 

T = V -\- _ 

'a ' oc ' 


AT, 


(AlO) 


Because the magnetic fields may be neglected, the forces 
FS = qaE2 contain only the electric field, their magni¬ 
tudes being unambiguously defined by the positions of par¬ 
ticles and boundary conditions. In the relativistic case, the 
scheme (lAlOII cannot be used directly, because F^ includes 
the Lorentz force gaV^ x depending on the velocity Vn 
determined at time r„ = nAT. However, in the leapfrog 
scheme, the particle velocities are defined at half-integral 
times Tn+i /2 = (n-|- l/2)Ar. The natural solution allowing 
us to retain the simplicity of the leapfrog scheme in this situa¬ 
tion is the application of Lagran ge’s interpolation formula for 


computing from the three values of the velocity (va 


%-\/2 


■^+l/2 


-rtT, + l 

Va = 


3g2-3/2 _ i5g;j+i/2 


(All) 


Thus, the complete integration scheme of the equations of 
motion takes the form: 


-#^+l/2 _ -#1—1/2 

Vet VoL 


+ xH")AT, 


.yn-ri /2 _ 


-# 1 + 1/2 




= - 


- lOUa 


—#1 + 1 —#1 I —I 

rot = 

,7i-3j2 _ ^pjwTi—1/2 


(A12) 


The fields and Fl^ acting on the particle are determined 
from the magnitudes of the mesh fields in the eight adjacent 
nodes f^ : 


is = ETki^ - wm - 0(1 - o 
+ ir+i,,',o(i-0(1-0 
+ iOi,;c(i-«'?K(i-0 
+ io+i(i-«^r)(i-o«'S 
+ ir+i,,'+i,o<(i-«^s) 

I rnTl /1 71 \ n 71 

-I- Eij+i^k+i(i - Wi )Wj Wk 

, ^71 71/1 71 \ 71 

+ ++l,j.fc+lWi (1 - Wj )Wk 

I rnTi n n n 

“t“ '^j '^k ■ 
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